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Abstract 



We review some well posed formulations of the evolution part of the Cauchy 
problem of General Relativity that we have recently obtained. We include also a 
new first order symmetric hyperbolic system based directly on the Riemann tensor 
and the full Bianchi identities. It has only physical characteristics and matter 
■ sources can be included. It is completely equivalent to our other system with these 

O ■ properties. 

O 

o\ '. 1 Introduction 

We will review some recently obtained well posed formulations of the evolution part of 
the Cauchy problem (see 0) in General Relativity, considered as the time history of the 
two fundamental forms of the geometry of a spacelike hypersurface, its metric g and its 
extrinsic curvature K. On such an hypersurface, for instance an "initial" one, these two 
quadratic forms must satisfy four equations, called constraints. The constraints can be 
posed and solved as an elliptic system by known methods. 

The proof of the existence of a causal evolution in local Sobolev spaces of g and K 
into an Einsteinian spacetime does not however result directly from the equations giving 
the time derivatives of g and K in terms of space derivatives of these quantities in a 
straightforward 3+1 decomposition of the Ricci tensor of the spacetime metric, which 
contains also the lapse and shift characterizing the time lines. These equations do not 
appear as a hyperbolic system for arbitrary lapse and shift, in spite of the fact that their 
characteristics are only the light cone and the time axis (see fll3"l). 



In this paper, we review different methods used recently for obtaining a hyperbolic 
evolution system for these geometrical unknowns which manifest their propagation gov- 
erned by the light cone. Particularly interesting are first order symmtric hyperbolic sys- 
tems, since they seem to be the most amenable to numerical computations. We construct 
explicitly such a system, using a harmonic time slicing condition. Its characteristics are 
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the physical light cone and the time direction orthogonal to the space sub-manifolds of 
the chosen slicing. Among the propagated quantities is, in effect, the Riemann curvature. 
The space coordinates and the shift are arbitrary and, in this sense, the system is gauge 
invariant. In fact, we can make the lapse also take arbitrary values by introducing a 
given arbitrary function in its definition. 



2 3+1 Decomposition of the Riemann and Ricci Ten- 
sors 

We suppose in all this paper that the spacetime is a smooth manifold, V = M x -ft, 
endowed with a metric g of signature (—,+,+,-1-), such that the "slices" are spacelike. 
This hypothesis is no restriction for globally hyperbolic spacetimes. 

We choose on V a moving coframe such that the dual frame has a time axis orthogonal 
to the slices M t while the space axes are tangent to them, namely we set 

9° = dt 
6 l = dx i + ftdt 

with t G 3? and x\ i = 1,2,3 local coordinates on M. The Pfaff derivatives d a with 
respect to 9 a are 

at 

In this coframe, the metric g reads 

ds 2 = g a(3 9 a 9? = -N 2 (9 ) 2 + ^ W. 

The t-dependent scalar iV and space vector (3 are called the lapse and shift of the slicing. 
Any spacetime tensor decomposes into sets of time dependent space tensors by projections 
on the tangent space or the normal to M t . 

We define for any t-dependent space tensor T another such tensor, doT, of the same 
type by setting 

A 9 r 

where Lp is the Lie derivative on M t with respect to f3. 

The extrinsic curvature of the slices is the time dependent space tensor given by 

K a = -\ N ~ l ®Q9ij 

We denote by an overbar geometric quantities associated with the Riemannian metric 
g = g,ijdx % dxi induced by g on the slices. 
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The Riemann curvature tensor of g admits the following decomposition into time 
dependent space tensors (cf. zero shift; 0, arbitrary shift; |19| , spacetime viewpoint) 



JV« = RAi ~ K> h Ku + KfK, 



ki 



ki ) 



From these formulas result the following ones for the Ricci curvature, where we have 
denoted by H the mean extrinsic curvature of the space slices, H — K\: 



r.. = -N- 1 d K ij +HK t] -2K im Kf -N^V^N + R ij 

Rot = N(-V h K^ + djH) 

R 00 = N(V%N -NK t] K lj +d H) 

From these identities results the following one for the Einstein tensor, leading to the 
so-called "Hamiltonian constraint": 

Goo = Roo ~ \gmR = \n\R - K.K + H 2 ) 

3 Wave Equations for K. Local Existence Theorems 

Using the relation between d^g^ and K^, we obtain (with no factor | in the index 
symmetrization) 



= -NV {i V h K j)h - V h d (i NK j)h - 2NR h ijm KJ? - NR^Kfi 
VhS h {NKij) + Vjdi(NH). 



We now use the expression for R 0i and Rij to obtain the identity 
= doRij — V(j-Rj)o 

= -^N-%K tj ) + d (HK tJ - 2K m Kf) - d^N^V faN) - 
NVidjH - V {l {K j)h O h N) - 2NR h ijm K™ - NR m{t K% + 
VhV^NKi^+HVjdiN 

This identity shows that for a solution of the Einstein equations 

the extrinsic curvature K satisfies a second order differential system which is quasi diag- 
onal with the principal part the wave operator, except for the terms VidjH. The other 
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unknowns g and N appear at second order except for the term c\)V jdiN '. These facts 
have been used in two different ways to obtain for K a quasi diagonal, wave type system. 

1. One eliminates at the same time the third derivatives of N and the second deriva- 
tives of H (cf. in the case of zero shift |J and for arbitrary shift |T(|) by requiring N to 
satisfy 

(N') d N + N 2 H = 

The second order equation for K reads then as a wave type nonlinear system: 

(K') NUK l3 = NQij + 9,, 

where we set 

nK .. = -N-^dodoKij + V^hKij 

NQij = -K lJ d H + 2g hm K m{l d K j)h + 4Ng hl g mk K lk K im K jh + 
(2V {i K j)l )d l N - 2HN- l diNdjN - 2d (i Nd j) H - 
WhN^Kij - K^VhN - N~ 1 K ij d h Nd h N + 
N- 1 K h{l d 3) NO h N + {V {i d h N)K j)h + 2NR h ljm K™ + 
NR m{i KY) - 2HV j 8 i N 
<£>ij = d pij — V(iPj)o 

The equation (JV') expresses that the time coordinate is harmonic. 
Using the expression for H we see that (N r ) reads 

d log(N/(detg) 1 *) = 

The general solution of this equation is 

N = or l {detg^) 

where a is an arbitrary tensor density such that 

d a = 0. 

The algebraic expression of the harmonic time-slicing condition is called "algebraic gauge." 

If we replace N by the value obtained above in the second order equation for K 
we obtain a quasi diagonal system with principal part the wave operator, with terms 
depending on g and its derivatives of order < 2. This system reduces to a third order 
hyperbolic system when we replace K by —(2N)~ 1 dog. The local existence theorem of 
Leray for the solution of hyperbolic systems gives immediately the local in time existence 
of solutions of this reduced system, in local in space Sobolev spaces, with the smallest 
index known for generic solutions of Einstein's equations, and domain of dependence 
determined by the light cone. It can be proved (see ||) that a solution of the reduced 
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system on M x / is a solution, in algebraic gauge, of the full Einstein equations if the 
initial data satisfy the Einstein equations on M . 



Remark 1. If we take for a an arbitrary given function we still obtain a quasi diagonal 
system for K, with additional a priori given terms. Such a generalization is equivalent 
to replacing (N f ) by the equation 

(N") N^doN + NH = f 

with / an arbitrary function. 

2. Replace in the term VjdiH the mean curvature H by an a priori given function h (a 



procedure used in jLj with h = 0, in the asymptotically euclidean case). The equations 
Qij = Qij become then, when N is known, a quasi diagonal second order system for K 
with principal part the wave operator, namely: 

(K") UK i5 = Pij + Gy 



where 



OKij = -dodoiN^Kij) + V h V h (NK i: 



Here, depends only on K and its first derivatives, on g, N and 8qN together with 
their space derivatives of order < 2. It is given by 



p.. = d (-HK ij + 2g hm K {im K^) + d (N- 1 V j d i N)+V h (d (i NK j)h ) + 
2NR h ijm K™ + NR m{i K% - HVfiN + NV&h, 



while Gij, zero in vacuum, is 

Gij = d pij — V(iPj)Q. 
When p, N and the sources p are known the above equation together with 
(<?') d 9ij = -2NKij 

are again a third order quasi diagonal system for g, hyperbolic if iV > and g is properly 
Riemannian. 

The condition H = h is a "mean curvature" gauge choice which imposes, through the 
equation Rq = that N satisfy the following elliptic equation on each slice 

(JV) V%N - (KijK ij - p° )N = -d h 

Note that for energy sources satisfying the strong energy condition we have — p° > 0, 
as well as K.K = K^K 1 ^ > 0, an important property for the solution of the elliptic 
equation. In the given mean curvature gauge the unknowns iV and g satisfy, for any 
choice of the shift /3, a mixed elliptic and hyperbolic system. 

Local existence theorems (i.e. in a neighborhood of M in M x dt) have been proven 



see [111]]) when M is compact and when M is asymptotically euclidean. 
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4 Case of Compact M 



Theorem 1. Let (M, e) be a smooth compact Riemannian manifold. Let there be given on 
M x I, with I = [0, T], a "pure space" smooth vector field (3 and a function h such that 

he f] C 3 - k (I,H k ), d h>0, d h^0. 

2<fc<3 

There exists an interval J = [0, £], £ < T such that the system {K"), (g'), (N) has one 
and only one solution on M x J 
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0<k<2 

with N > and g uniformly equivalent to e, taking the initial data 

9ij(0,.) = <y tj E H 3 
#y(0,.) = k l3 EH 2 

if 7 is a properly riemannian metric uniformly equivalent to e and k.k ^ 0. 

5 Case of (M,e) Euclidean at Infinity 

The manifold M is the union of a compact set and a finite number of disjoint sets (its 
"ends") diffeomorphic to the exterior of a ball in 3? 3 ; the smooth given metric e reduces 
on each end to the euclidean metric. The weighted Sobolev space H s $ on (M, e) is the 
completion of in the norm 

\\f\\H., t = (f Z 0<k< y k+25 \D k f\^(e))* 
jm — 

where a 2 = 1 + d 2 and d is the distance in the metric e to some fixed point in M . 

The use of the theory of elliptic equations on an asymptotically euclidean manifold 
in H s s spaces (see 0) and of weighted energy estimates leads to the following theorem: 

Theorem 2. The system (1),(2),(3) with Cauchy data 7 and k on the manifold (M, e) 
euclidean at infinity with given j3 and h on M x [0, T] has one and only one solution 
(g,K,N) on M x J, J = [0,£] a sufficiently small subinterval of I, if the Cauchy data 
are such that 7 — e G #3,-1 and is uniformly equivalent to e while k G H 2t o and k G H\ t \. 
The solution belongs on each M t to the same functional spaces as the data. It is such 
that N > and g is uniformly equivalent to 7. 
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Remark 2. The estimates show in fact that g — 7 is such that for each t we have: 

g t - 7 G #3,0 

which is a stronger asymptotic fall off than 7 — e or g t — e: this property is related to the 
Arnowitt-Deser-Misner theorem of mass conservation. 

6 First Order System (Vacuum) 

We use the wave equation satisfied by K to show (see flQfl , that in spacetime di- 
mension 4 a solution of the vacuum Einstein equations, together with the harmonic time 
gauge condition, satisfies a first order symmetric system, hyperbolic if g is properly Rie- 
mannian and N > 0. Such a system could be useful to establish a priori estimates relevant 
to global problems. It may be important for numerical computations (see ||) because 
symmetric hyperbolic systems occur in many areas of mathematical physics, in particular 
in fluid dynamics, and effective codes have been developed to study such systems. 
We have obtained for the unknowns g, K, N the equations 

(1) d o9ij = —2NKij 

(2) d N = —N 2 H 

(3) UKij = Q i:j 

To obtain a first order system we take as additional unknowns: 

N^doKij = Uj 
V h K i3 = M hij 
dilogN = cii 
N^dodilogN = a 0i 
N^VjdilogN = ctji 

We take as equation (3') 

(3') doK i;i 
The equation (3) gives 

(4) doLij 
In three space dimensions the Riemann tensor is a linear function of the Ricci tensor: 

■Rlijm = QljRirn Qim-Rjl QijRlm 9lm-Rij "^\9lj9iTn 9ij9lrn)R 

Using the equation Rij = to express Rij, we write Qij as a polynomial in the unknowns 
and g % K 



NL; 



NV h M hi] = N(HL { j - Q t 
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By using the identity 

doF % = V h (NKij) - V {i (NK*) 
we see that for an arbitrary covariant vector Ui we have 

doV h Ui = V h d oUl + Ul (V h (NKl) - Vi(NK l h ) - V l {NK ih )) 
From an analogous formula for tensors we deduce the equation 

(5) d M hij - NV h Lij = N[a h L l3 + (M^K^ + K^M^ - K l{j M l i)h )+ 

Ki^K^au + Oi)K l h - a l K i)h ) 

On the other hand (2) implies 

(6) Boat = Na 0l = -N(M lh h + H<n) 
while dhi and must satisfy 

(7) d a hi - NV h a 0l = JV(a,(M w ' + M ih l - M l ih ) + ai{K\a h + K l h a t - a l K lh ) + a h a 0i ) 

The condition (2) together with the Einstein equation imply that N satisfies the 
inhomogeneous wave equation (this step brings in the Hamiltonian constraint determined 
by Goo) 

d d N - N 2 V h V h N = —N 3 KijK ij + 2N 3 H 2 

Hence by differentiation, use of the Ricci formula, the definitions of au and ao« and 
simplification through the use of R% = 

(8) d a 0l - NV h a hi = N(2a h a ih - 2K hl M i hl + HM ih h - R>ta h ) + 

N ai (H 2 + 2a h h + 2a h a h - 2K hj K h i) 

We use again the equation Rij = to replace R\j by its value in terms of the un- 
knowns. We have obtained a first order system (1,2, 3', 4, 5, 6, 7,8) in all the unknowns. The 
right hand sides are polynomial in the unknowns and they do not depend on their 
derivatives. The left hand sides are linear operators on all the unknowns. Their coeffi- 
cients depend on these unknowns, and not on their derivatives except for the derivatives 
of g which appear through the Christoffel symbols T^. We can use the identity given 
above to write their evolution equation 

d Q Y% = N(M% - M (y 5 - a h K l3 - a^). 

However, if we wish to display an explicitly covariant quasi linear first order system for 
all the unknowns, in particular if M is not diffeomorphic to 3? 3 , we introduce on M an 
a priori given connection E, possibly derived from a metric e which will also be used to 
define the Sobolev spaces on M. The derivatives in the connection of g and the given 
connection differ by terms linear in a tensor S, the difference of the two connections. The 
evolution of this tensor is given by 

d S* = N(M% - M {i .* - a h K i3 - a {i K*) - d E*. 
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We have now obtained a quasi linear first order system covariant for all the unknowns. 
Its characteristic matrix, obtained by replacing in the principal matrix the operator d 
by a covariant vector £, consists of blocks around the diagonal, some reduced to one 
element £o and some 4x4 matrices with determinant £g — iV 2 £ l £j. The characteristics 
are the light cone and the axis orthogonal to the time slices. On the other hand, the 
system can be symmetrized by multiplication by a matrix consisting of blocks around 
the diagonal equal to one element, 1, or the matrix (g^). In other words, the system is 
a first order symmetrizable hyperbolic system, with domain of dependence determined by 
the light cone. Known local existence theorems apply to such a system. 



7 Bianchi Equations 

Instead of taking as new unknowns the first derivatives of K to obtain a first order 
symmetric hyperbolic system, one can use the Riemann tensor of the spacetime metric, 
which is linear in these derivatives, to rewrite the system in Section 6. It satisfies the 
Bianchi identities: 

These identities imply by contraction and use of the symmetries of the Riemann tensor 

VaR^fa + V 7 i?/3 A1 + VpRjfj, = 0. 

If the Ricci tensor R a p satisfies the Einstein equations 

Ra/3 — Pa/3 

then the previous identities imply the equations 

The first equations with (a/37) = an d the last one with /x = do not contain 

derivatives of the Riemann tensor transversal to M t ; they are considered as constraints. 

In a particular case (see p0| ) of the general system considered below, the conservation 
of an initially vanishing Riemann tensor in vacuum was obtained. That result does not 
require analyticity when it is written in symmetric hyperbolic form as below. 

We wish first to show that the remaining equations are, for n = 3 in the vacuum 
case, when g is given, a symmetric first order hyperbolic system for the double two-form 
R a(3,\n- For this purpose, following Bel (see ||) we introduce two pairs of "electric" and 
"magnetic" space tensors associated with a spacetime double two-form A, 

N 2 Eij = A 0i fij 
Dij = —ilihkVjimA hk,lm 

NH l3 = ^Vih k A hk fij 

NB 3l = \r) lhk A^ hk 
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where 77^ is the volume form of g. It results from the symmetry of the Riemann tensor R 
with respect to its first and second pairs of indices (R is a "symmetric double two-form" ) 
that if A = R, then E and D are symmetric while = B^. The Lanczos identity (see 



15 1) for a symmetric double two form like R, with a tilde representing the spacetime 



double dual, is 

where C Q/ g = R a p — \g a pR. It follows that when R Q/3 = \g a p, then E = —D and H = B. 
In order to avoid introducing unphysical characteristics, and to be able to extend the 
treatment to the non vacuum case, we do not use these properties in the evolution 
equations, but write them as a first order system for an arbitrary double 2-form A as 
follows 

Vo^hfc,oj + VfcAo^oj — ^hAok,oj = 
V o A o i)0j + V h A\ Qj = V oPj i - Vp 0l 

and analogous equations with the pair (Oj) replaced by (lm). One obtains a first order 
system for the unknowns E, H, D, B by using the following relations deduced from the 
definition of these tensors 

Ahk,oj = Nrf hk Hij 

Ahk,lm = hkrflmPii 

A^^hk = Nrf hk Bji 

The system obtained has a principal matrix consisting of 6 identical 6 by 6 blocks around 
the diagonal, which are symmetrizable and hyperbolic. Hence, the system is symmetric 
hyperbolic, when g is a given metric such that g is properly Riemannian and iV > 0. 
To relate the Riemann tensor to the metric g we use the definition 

doQij = -2NK i:j 

We consider two possibilities: each introduces a choice of gauge for N. 
1. We fix the mean curvature of the space slices, i.e. we set 

H = h 

with h being a given function. We deduce from the identity giving R t ° kl the "elliptic" 
identity 

V h V h K t , - R)K hi - R hljm K hm = V h (NR\ jh ) - V^JV" 1 ^) + V^H. 

We write an elliptic system linking the symmetric tensor K with a double 2-form A when 
H = h by symmetrization of the above identity, replacement of R by A, H by h, and the 
spacetime Ricci tensor by p, zero in vacuum. Namely, we consider the system 

V h V h K t3 - fyK j)h - R hijm K hm = V h {NA\ i]j)h ) - V {l {N- l Pj)o ) + V&h. 
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This system is to be solved globally on each slice M t , like the equation for N in mean 
curvature gauge. We have obtained a mixed elliptic-hyperbolic system for the unknowns 
g, N, K, E, D, H, B. It can be proved that this system is well posed in the case of compact 
or asymptotically euclidean M, as long as the elliptic operators are injective. This prop- 
erty holds for metrics g in a neighborhood of a metric 7 for which it holds, for example 
if M = 9ft 3 and 7 is the euclidean metric, or if M = S 3 and 7 is the canonical metric of 
S 3 . 

2. We choose the algebraic form of the harmonic time-slicing condition, that is we set 
N = a- l \detg\ 1 ' 2 , with a a given tensor density. We follow the ideas used by Friedrich 
(see ||13|| ) for the Weyl tensor to write a symmetric hyperbolic system for K and F, 
namely we consider the following identities, deduced from the definition of K and the 
3 + 1 decomposition of the Riemann tensor: 

+ NV h K l0 = K tJ d h N - KfaN - R 0( J 

and 

doKij + NRij + VjdiN = -2NR\ 0j - NHK tJ + NR iy 

We obtain equations relating F and K, for a given double 2-form A, by replacing in these 
identities the Riemann tensor by A and the Ricci tensor of spacetime by a given tensor p, 
zero in vacuum. To deduce from this system a symmetric hyperbolic first order system, 
with the algebraic form of the harmonic gauge, one uses the fact that in this gauge one 
has 

T h lh = dilogN + diloga 
The second set of identities leads then to the following equations: 

doK.j+Nd^ = N{V2t) m -{t^ h -d i ^){V l jl -d j loga)) + 

NVjdiloga - 2NA\ 0j - NHK i5 + N Plj 

The system obtained for K and F has a characteristic matrix composed of 6 blocks 
around the diagonal, each block a 4 by 4 matrix that is symmetrizable hyperbolic, with 
characteristic polynomial ^ a ^a- 

The whole system for A, K, F, g is symmetrizable hyperbolic, with characteristics the 
light cone and the normal to M t . 

Remark. It is somewhat involved to prove that a solution of the constructed system 
satisfies the Einstein equations if the initial data, 7, k, satisfy the constraints, but we 
can argue as follows. We consider the vacuum case with initial data 7 and k satisfying 
the Einstein constraints. These initial data determine the initial values of F, and also, if 
f3 and iV are known, the initial values of Aij thm , Aj h ^ , A i0 j h by using the decomposition 
formulas. (We set A equal to the Riemann tensor on the initial surface.) We use the 
Lanczos formula to determine A iQt j Q initially. We know that our symmetrizable hyperbolic 
system has one and only one solution. Since a solution of Einstein's equations in algebraic 
gauge, proved to exist in a preceding section, satisfies together with its Riemann tensor 
the present system and takes the same initial values, that solution coincides with the 
solution of the present system in their common domain of existence. 
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8 A Non-Strict Hyperbolic System with Arbitrary 
Lapse and Shift 

Lemma. The following combination of derivatives of components of the Ricci tensor of 
an arbitrary spacetime is a quasi diagonal hyperbolic operator for the extrinsic curvature 
of the space slices: 

= d d Rij — <9 V(j-Rj) + VjdiR 00 

Proof. A straightforward computation using previous results shows that 

Ajj = d nK ij + dod (HK ij -2K im KJ l )-d d (N- 1 V j d i N) + 

d (-V (i (K j)h d h N) - 2NR h ijm K™ - NR m{i K^ + HV &N) + 
Vjdi(NAN - N 2 K.K) + C i5 

with 

aK .. = -d^N-'doK^) + V h V h (NK, tj ), A = V,V A 

and 

dj = VMNdoH) - doiNVjdiH) 

We see that Cij contains terms of at most second order in K (and also in N) and first 
order in g (replace d gij by —2NKij). 

The identity given above shows that for a solution of the Einstein equations the 
extrinsic curvature K satisfies, for any choice of lapse and shift, a third order differential 
system which is quasi diagonal with principal part the hyperbolic operator BqO. The 
other unknown g appears at second order except for terms appearing through VjdiAN. 

The system for g and K is not hyperbolic in the usual sense of Leray because of the 
third order space derivatives of g appearing in Vj-c^AiV. 

Theorem 3. The system for g, K is equivalent to a system hyperbolic non strict in 
the sense of Leray-Ohya with local existence of solutions in Gevrey classes and domain 
of dependence determined by the light cone. 

Proof. Replace in the equation K by — (2A) -1 <9o(?: this gives a quasi diagonal system 
for g, but with multiplicity 2 for do in the principal operator. 

The system for g, K can be turned into a hyperbolic system by a gauge choice (or 
redefinition of N) leading either to an elliptic or to a hyperbolic equation for N containing 
an arbitrary function. Here, we sketch a forthcoming more detailed treatment (see 0). 

1. We give arbitrarily a constant c, taken to be zero in the asymptotically flat case 
and strictly positive if M is compact. We choose arbitrarily a smooth function / on 
M x 3?, with appropriate asymptotic behavior if M is euclidean at infinity and / < and 
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/ ^ on each M t if they are compact. We require N to satisfy on each M t the elliptic 
equation depending on g 

AN-cN = f 

and we determine No = doN and N 0Q = <9q A" from the elliptic equations obtained by 
differentiating the equation for N and replacing dog by —2NK. It can be shown that 
the mixed hyperbolic and elliptic system obtained for g, K, N, N , N 00 is well posed, in 
appropriate functional spaces. 

2. We choose arbitrarily a smooth function / on M x dt and require N to satisfy the 
wave equation. 

N~ 2 d d N -AN = f 

We use this wave equation to reduce the term d^N" 1 VjdiN) in A^-. It can be shown that 
the system for g,K,N, N is hyperbolic, hence well posed in local Sobolev spaces, if N 
is determined by the hyperbolic equation obtained by differentiating the wave equation 
satisfied by N. 
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